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Abstract 

A universe consisting of two interacting perfect fluids with the same 
4-velocity is considered. A heuristic mean free time argument is used 
to show that the system as a whole cannot be perfect as well but nec- 
cessarily implies a nonvanishing bulk viscosity. A new formula for the 
latter is derived and compared with corresponding results of radiative 
hydrodynamics. 
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1 Introduction 



In the realm of cosmology bulk viscosity is the most favorite dissipative phe- 
nomenon. Different from shear viscosity and heat conductivity, it is com- 
patible with the symmetry requirements of the homogeneous and isotropic 
Priedmann Lemaitre Robertson Walker (FLRW) universes. In the simplest 
cosmological models there is no way to study entropy producing processes 
except throTigh bulk viscosity. 

While this is obvious on formal grounds it is probably fair to say that 
the degree of 'understanding' bulk viscosity physically is by far less than 
that for shear viscosity or heat flux. Although the corresponding effect for a 
simple gas is known since the work of Israel (1963) and a radiative bulk vis- 
cosity coefficient was derived by Weinberg (1971) who also gave an analysis 
on its role in cosmology, followed by subsequent investigations of Straumann 
(1976) and Schweizer (1982), there remained the desire of some more intu- 
itive insight beyond the involved calculations of radiative hydrodynamics. 
This problem was addressed in the last part of a paper by Udey & Israel 
(1982) who argued that for a two-fluid universe the mechanism responsible 
for bulk viscosity is a microscopic heat flux that compensates the temper- 
ature differences caused by different cooling rates of the two components. 
Their consideration of this point was based on the following semiquantitative 
argument. Let r be the characteristic time for the interaction between both 
fluids such that during a time interval r the perfect fluid components may 
be considered as effectively insulated from each other, resulting in different 
adiabatic cooling rates due to their different equations of state. 

The present paper was inspired by this kind of arguing. While we shall 
confirm below that the difference in the cooling rates is indeed the essential 
point, we shall avoid in our investigation the introduction of microscopic 
heat fluxes. Although this concept may be helpful on small scales, it seems 
less convincing on large scales, e.g., of the Hubble scale. The microscopic 
gradients at different points had to conspire in order to produce a nonvanish- 
ing bulk viscosity in homogeneous and isotropic universes, which, however, 
is incompatible with the symmetry requirements of the latter. Moreover, the 
apparent reduction of a bulk viscous pressure to heat fluxes is not consistent 
with the fact that both phenomena arc basically independent. 

It is the aim of this paper is to show that different cooling rates for two 
perfect fluids are sufficient for the existence of a nonvanishing bulk viscosity 
of the system as a whole. No additional concept like that of a heat flux over 
intermolecular distances has to be used. The basic idea is to study a universe 
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of two different interacting perfect fluids and to ask for the conditions under 
which an effective one-fluid description is possible. It turns out that this 
one-fluid universe is neccessarily dissipative. 

The present paper is semiquantitative throughout. We are not claiming 
to improve any of the technically rather complicated calculations in radiative 
hydrodynamics. Our objective is to achieve a kind of phenomenological 
'understanding' of the bulk pressure phenomenon in the expanding Universe. 
On this level of description, we shall derive a new formula for the coefficient 
of bulk viscosity in a two-fluid system. 

In section 2 the relevant relations for two noninteracting perfect fluids 
are presented and the corresponding cooling rates in an expanding universe 
are obtained. Section 3 is devoted to an effective one- fluid description for 
two fluids with mutual interaction. An explicit expression for the coefficient 
of bulk viscosity is derived with the help of a mean free time argument. 
The latter result is compared with work on radiative hydrodynamics in the 
context of relativistic kinetic theory in section 4. Section 5 gives a brief 
summary of the paper. 



2 Two-fluid dynamics 

The content of the Universe is assumed to be describable by an energy mo- 
mentum tensor T*^ that is the sum of two different perfect fluid contributions 
which share the same 4-velocity: 

with (A = 1, 2) 

T^ = PAu'u^+PAh''' . (2) 

PA is the energy density and pA is the equilibrium pressure of species A. 
is the common 4-velocity and /i*^ is the projection tensor h^^ = g"^^ + n'u^. 
Let us first deal with the case that the energy momentum conservation laws 
hold for each fluid separately: 

T%k = ^. (3) 

implying the energy balances 

PA = -Q{pa+Pa) , (4) 



3 



with the fluid expansion = n*j. pA = PA,iU^ etc. Because of (^) both 
fluids evolve independently except for their mutal gravitational coupling. 
The particle flow vector N\ of species A is defined as 

N\ = riAu' , (5) 

where ua is the particle number density. Particle number conservation is 
expressed by N^.^ = 0, equivalent to 

ua + QnA = . (6) 

Let us further assume equations of state in the general form 

PA=PA{nA,TA) (7) 

and 

PA = PA {nA,TA) , (8) 

i.e., let particle number densities ua and temperatures Ta be our basic 
thermodynamical variables. The temperatures of both components will be 
different in general. 

Differentiating relation (^), using the balances (Q) and (^) as well as the 
general relation 

dpA ^ PA+PA Ta dp a 

duA riA uaQTa 
that follows from the requirement that the entropy is a state function, we 
find the following expression for the temperature behaviour: 

T. = -r.e«i . (10) 

dpA/oTA 

It is obvious that the temperatures of both components behave differently 
for different equations of state. With Q = 3R/R, where R is the scale 
factor of the Robertson- Walker metric, the equations of state pi = riikTi, 
pi = SnikTi reproduce the well known Ti ~ R~^ behaviour for radiation. 
With p2 = n2kT2, p2 = n2mc^ + ^n2kT2 one obtains T2 ~ R~^ for matter. 

3 Effective one-fluid dynamics 

The hitherto independent fluids are now allowed to interact. We try to find 
an effective one fluid description for the Universe as a whole, characterized 
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by the particle number density n = ni+n2 and an equilibrium temperature 
T. The overall equations of state are 



p = p{n,T) (11) 

and 

p = p{n,T) , (12) 

where p is the equilibrium pressure and p is the energy density of the system 
as a whole. The equilibrium temperature T is defined by (cf. Udey & Israel 
1982) 

pi{ni,Ti) + p2{n2,T2) = p{n,T) . (13) 
As we shall show below this implies 

pi{ni,Ti) + p2{n2,T2) p{n,T) . (14) 

For perfect fluids the difference between both sides of the latter inequality 
is the viscous pressure vr Q 

IT = pi{ni,Ti) + p2{n2,T2) - p{n,T) . (15) 

The existence of a nonvanishing viscous pressure is a consequence of the 
different temperature evolution laws of the subsystems. This is most easily 
understood by the following simple mean free time argument. Let r be the 
characteristic mean free time for the interaction between both components. 
The time r is assumed to be much larger than the characteristic interac- 
tion times within each of the components. Consequently, the latter may 
be regarded as perfect fluids on time scales of the order of r. The inter- 
action between the fluids is modelled by 'collisional' events, where r plays 
the role of a mean free time between these 'collisions'. During the time 
interval r, i.e., between subsequent interfluid interaction events, both com- 
ponents then evolve according their internal perfect fluid dynamics, given 
W (@)) ® s-iicl (10). Assume that through this interaction an element of 



the cosmic fluid is in equilibrium at a proper time rjQ at a temperature 
T im) = Ti (rjo) = T2 (rjo) with p (%) = pi (%) + P2 (%)• Here, p (%) and 
PA ivo) are shorts for p [n (%) , T (r/o)] and pA [riA (%) , Ta (??o)], respectively. 
Using the condition (|l3|) at the proper time r/o + r up to first order in r, i.e.. 



PAivo + r) = pA{vo) + TPAivo) + ■■■ (16) 



^For real fluids this is not neccessarily true (see, e.g. Kirkwood & Oppenheim 1961 
chapter 7.7) 
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and 

P iVo + r) = p{i]o) + TP (rio) + ... (17) 

where pA (%) = Pa [nA (%) , Ta (%)] and p (770) = p [n (r/o) , T (%)] , applying 
(^) on the l.h.s. of ( p^ ) and the general relation 

dp p + p T dp 



dn n n dT 



(18) 



on its r.h.s., one finds 



at the point r] = rjQ. 

During the following time interval r, i.e., until a subsequent 'collision', 
the subsystems move freely according to their proper dynamics given by 
(I), (I) and (|10|). At the time rjo + t we have T (r?o + r) / Ti (r/o + r) / 
?2 (% + ''") 7^ ^ (^0 + ''") in general. With 

Ta{vo + t) =TA{vo)+TfA{vo) + ... (20) 

and 

r(r?o + T) =r(r?o)+Tr(r?o) + - (21) 
we find from (|l^) and ([K 



dpi/dT 


dp2/dT^ 


dpi/dT 


dp2/dT, 


dpi/dT 


dp/dT\ 


dpi/dT 


dp/dT ) 


dp2/dT 


dp/dT\ 



T,-T=-reTf*l^-*^l , (23) 



up to first order in r. Due to the different cooling rates ( |To| ) and (1£) there 
occur temperature differences at any point of the expanding fluid: Both 
differences in the temperatures of the components, i.e., between Ti and T2, 
and differences between the temperature of each of the components and the 
temperature T of the system as a whole. 

In order to arrive at our conclusion vr 7^ one has simply to consider the 
sum of the partial pressures pi and p2 at = 7/0 + t up to first order in r. 
The latter may be written as 

pi{ni,Ti) +p2{n2,T2) = pi (ni, T) + p2 (^^2, T) 

dpi dp2 



+ {Ti-T)^ + {T2-T)-^ . (25) 
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Inserting the temperature differences Ti — T and T2 — T at tjq + t from (|2 
and dH), we find 

dpi fdpi/dT dp/dT\ 



pi (ni,Ti) +P2 ("-2,72) = p{n,T)-TTQ 



dT \dpi/dT dp/dTj 
dp2 (dp2/dT dp/dT\ 
^ W \dp2/dT ~ dp/dTj 



(26) 



where p (n, T) = pi (ni, T) +p2 (n2, T) was used. Applying the zeroth-order 
relations dp/dT = dpi/dT + dp2/dT and Sp/ST = dpi/dT + dp2/dT in 
the bracket on the r.h.s. of ( p^ ) one gets, after a simple rearrangement, the 
first-order result 

pi {ni,Ti) +P2 (n2,T2) =p{n,T) + TT , (27) 

where 

/Opi/ar _ dp/dT \ f dp2/dT _ dp/dT \ 

^ ^ dT\dpi/dT dp/dT J \dp2/dT dp/dT J ^ ' 

is generally different from zero. This proves our initial statement that a 
system of two interacting perfect fluids is not perfect as well. While the 
energy momentum tensors of the subsystems are given by (|2|), the system 
as a whole is characterized by 

T'' = pu^u^ + (p + vr) /i^^ . (29) 

To separate bulk viscosity from any other dissipative phenomenon we have 
ignored here the possibility of nonvanishing heat fluxes and shear stresses in 
inhomogeneous and anisotropic cosmological models. 
Prom the definition 

vr = -CG , (30) 
of the bulk viscosity Q, the latter is found to be given by 

i = -rT%{^^-^A{^^-^A , (31) 



where 

dpA _ f dpA\ _ dpA/dT dp _ f dp\ _ dp/dT 



dp A ydpA/nA dpA/dT ' dp \dpj ^ dp/dT 



(32) 
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For 'ordinary' matter dpA/dpA lies in the range 1/3 < dpA/dpA < 2/3 
(see the equations of state below ([lO|)). The lower limit corresponds to 
radiation, the upper one to matter, dp/ dp wih take a value intermedi- 
ate between dpi /dpi and dp2/dp2- Assuming without loss of generality 
dp/dp > dpi/dpi, we shall have dp/dp < dp2/dp2- Consequently, C ^ 0, 
i.e., the entropy production is positive which agrees with the second law of 
thermodynamics (see, e.g., de Groot, van Leeuwen & van Weert 1980). 

To the best of our knowledge the expression ( |3l| ) for the coefficient of 
bulk viscosity is a new result . Although based on heuristic arguments its 
structure is rather general. Formula ( |3l| ) is valid for general equations of 
state (0), (^) and ([iTl), (|T^). Until now we did not specify to the case that 
one of the components obeys the equations of state for radiation. 

4 Relation to radiative hydrodynamics 

The expression ( ^l|) for the coefficient of bulk viscosity is similar but not 
identical to the corresponding expressions of radiative hydrodynamics found 
by Weinberg (1971), Straumann (1976), Schweizer (1982), Udey &: Israel 
(1982) and Pavon, Jou & Casas- Vazquez (1983). The reason for this is a 
physical one and not due to the semiquantitative nature of our considera- 
tions. Both components of the system are treated as fluids with different 
equations of state in our setting. The final result (^) is therefore symmet- 
ric under a change of the labels 1 and 2 that identify both components. In 
the work of Weinberg (1971), Straumann (1976), Schweizer (1982), Udey & 
Israel (1982) and Pavon, Jou & Casas- Vazquez (1983) on the other hand, 
both components are treated asymmetrically. While one of the components 
is a fluid as well, the second one, a radiation component, is described with 
the help of kinetic theory. The main asymmetry lies in the assumption of the 
mentioned authors that the radiation component is allowed to deviate from 
local equilibrium while the fluid component is not. Of course, the result for 
the coefficient of bulk viscosity is not symmetric in the components either 
and a coincidence with ( ^l]) cannot be expected. 

It might be useful to compare some of the basic relations of radiative 
hydrodynamics with our framework. For definitenes, let component 1 of 
our analysis be the radiation and component 2 the material component. 
Since there does not appear a separate radiation temperature in the papers 
by Weinberg (1971), Straumann (1976), Schweizer (1982), Udey &: Israel 
(1982) and Pavon, Jou & Casas- Vazquez (1983), let us eliminate the latter 
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in (|T3D. With pi (ni,ri) = pi (ni,r2) + (Ti - T2) dpi/dT + ... we find, up 
to first order in the temperature difference, 



p {n, T) = pi (ni, Ts) + p2 (ng, T2) + p , (33) 

where 

n= (T, - To^ 

dT 



P^(Ti-T2)^. (34) 



Using (p^) with (p^) we have 



dT \dpi dp2 

For comparison with the results of Udey & Israel (1982) for radiative hy- 
drodynamics it is helpful to rewrite the latter expression as 

pui + tTQ— — , (36) 



^ ' dT \dp2 dp 

where 



9pi (dpi dp\ 
PUI = -tTQ— (^^ - ^ j • (37) 



With (|36| ) and (|37| ) relation ( |33D may be compared with formula (31) of 
Udey & Israel (1982), which, in our notation, reads 

p{n,T) = p{n,T2)+aT^B . (38) 

The quantity p in our equation (^) is the counterpart of the term aT^B 
in Udey Sz Israel (1982), which describes the deviation of the radiation 
component from equilibrium. With Tdpi/dT = AaT^ and dpi /dpi = 1/3 
for radiation and specification of (35b), (45) and (47) in Udey & Israel (1982) 
to the Eckart case (i.e., neglecting relaxation and cross effects), we find that 
aT^B in ( ^8| ) (equation (31) in Udey & Israel 1982) coincides with pm of 
our equation (|37|). The circumstance that there exists a difference between 
p and pui reflects the above mentioned fact that, different from radiative 
hydrodynamics, in our setting both components are allowed to deviate from 
equilibrium. For dp2/dp2 ~ dp/ dp this difference becomes negligible and 
we have p ~ pui. 

A similar statement holds for the difference of the temperatures T2 — T. 
By virtue of the zeroth-order identity 

(39) 



dp , 


/ dp 


dp2\ 


dpi 


f dpi 


dp2\ 


dT ' 


[d^~ 


9/02/ 


dT 


[dpi 


dp2) 
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the temperature difference (24) may be written as 

dpi/dT f dpi dp2 



^ ^ ^ dp/dT \dpi dp2 



where 



(41) 



dp/dT \dpi dp, 

if specified to radiation, coincides with Weinberg's relation (2.38) (Weinberg 
1971) for radiative hydrodynamics. The condition for T — T2 ~ (T — T2)^r 
is again dp2/dp2 ~ dp/ dp. 

Using the corresponding intermediate steps that led to (33) now for the 
sum of the pressures, we get the first-order relation 

pi (ni,ri) +P2 (n2,T2) = pi (ni,r2) + p2 (n2,T2) +p , (42) 

with 

P=m-r2)— = — p. (43) 

For the first two terms on the r.h.s. of (^) we may write, again up to first 
order, 

Pi{nuT2)+P2{n2,T2)=p{n,T2) = p {n,T) + {T2 - T) ^ . (44) 
From (^) and ( |35D one has 

T^-T^-{^y%. (45, 

and, consequently, 

pi (nijTi) +P2 (n2,?2) = p{n,T) + IT , (46) 

with 

/ dpi dp\ ^ 

This is relation (33) of Udey & Israel (1982). 

In order to avoid misunderstandings we point out that it was not nec- 
cessary for our derivation of formula (^) to introduce the quantities p and 
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p. These quantities are useful, however, for the comparison with work done 
in radiative hydrodynamics. 

By virtue of the identity p9|), from (|3l]) may be written as 



While it is generally not to be expected that this expression coincides with 
Weinberg's coefficient (Weinberg 1971) 

f„ = .T^l - P\' , (49) 



dT \dpi dp J 

it is obvious that both results are the closer the better the approximation of 
dp/dp by dp2/dp2 in ( |3T[ ) will be. In other words, the results (48) and ( |49| ) 



are similar if the matter component dominates the behaviour of the system 
as a whole. 

The bulk viscosity coefficients C and are related by 



. dpi 
dT 

Using (31) on the r.h.s of the latter equation we obtain the following relation 



/ dp 


dp2\ 


f dpi 


dp\ 


\d~p~ 


dp2) 


[dpi 





between C and Cw- 

^C = Cw, (51) 

where 

dp2 ^ dp2/dT 

dp dp/dT ■ ^ ' 

Relation ( pD shows again that the generally different expressions for and 
(^w will become similar if the overall behaviour of the system is determined 
by the fluid component 2. 

For the specific case of a mixture of radiation and matter with equations 



of state pi = nikTi, pi = SnikTi, p2 = n2kT2, p2 = n2mc^ + ^n2kT2 the 



bulk viscosity coefficient (|3^) reduces to 



C = ^mkT-^^ . (53) 
3 2ni + n2 

Since dp2/dp = n2/ (2ni + n2) one has C ~ Cw for ni <C n2- The lower 
the ratio ni/n2 of the photon number density to the number density of the 
matter particles the closer to unity is the ratio Cw/C- 
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A fluid description of tlie Universe makes sense as long as r <^ , 
where = 0/3 is the Hubble parameter. Combining ( |5^ ) with (|30| ) we find 
that this condition is consistent with | vr pa, i-e-, the magnitude of the 
nonequilibrium part tt of the pressure is much smaller than the equilibrium 
pressures as it is neccessary for a first-order approach like that of the present 
paper to be valid. 

5 Summary 

This paper is a heuristic attempt to clarify the origin of bulk viscosity in 
the expanding universe. Characterizing the interaction between two differ- 
ent fluids, each of them perfect on its own, by an effective mean free time 
parameter r and assuming a free evolution of both components according to 
their internal perfect fluid dynamics during the time interval r, i.e., between 
subsequent interfluid interaction events, the different cooling rates (due to 
different equations of state) of the components lead to a nonvanishing bulk 
pressure of the system as a whole. A new formula for the coefficient of bulk 
viscosity of a two-fluid mixture was obtained and the relation of this expres- 
sion to the results of radiative hydrodynamics was clarified. 
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